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Abstract 

We present theoretical predictions for the transverse-momentum distribution of 
J/ij: mesons promptly produced in 77 collisions within the factorization formalism 
of nonrelativistic quantum chromodynamics, including the contributions from both 
direct and resolved photons, and we perform a conservative error analysis. The 
fraction oi J/tp mesons from decays of bottom-flavoured hadrons is estimated to 
be negligibly small. New data taken by the DELPHI Collaboration at LEP2 nicely 
confirm these predictions, while they disfavour those obtained within the traditional 
colour-singlet model. 

PACS numbers: 12.38.Bx, 13.65.-|-i, 14.40.Gx 




Since its discovery in 1974, the J/'tp meson has provided a useful laboratory for quanti¬ 
tative tests of quantum chromodynamics (QCD) and, in particular, of the interplay of per¬ 
turbative and nonperturbative phenomena. The factorization formalism of nonrelativistic 
QCD (NRQCD) |l[] provides a rigorous theoretical framework for the description of heavy- 
quarkonium production and decay. This formalism implies a separation of short-distance 
coefficients, which can be calculated perturbatively as expansions in the strong-coupling 
constant from long-distance matrix elements (MEs), which must be extracted from 
experiment. The relative importance of the latter can be estimated by means of velocity 
scaling rules, he., the MEs are predicted to scale with a dehnite power of the heavy-quark 
{Q) velocity v in the limit u -C 1. In this way, the theoretical predictions are organized 
as double expansions in and v. A crucial feature of this formalism is that it takes into 
account the complete structure of the QQ Fock space, which is spanned by the states 
n = dehnite spin S', orbital angular momentum L, total angular momentum 

J, and color multiplicity c = 1, 8 . In particular, this formalism predicts the existence of 
color-octet (CO) processes in nature. This means that QQ pairs are produced at short 
distances in CO states and subsequently evolve into physical, color-singlet (CS) quarkonia 
by the nonperturbative emission of soft gluons. In the limit u —> 0, the traditional CS 
model (CSM) is recovered. The greatest triumph of this formalism was that it was 

able to correctly describe the cross section of inclusive charmonium hadroproduction 
measured in pp collisions at the Fermilab Tevatron 0 , which had turned out to be more 
than one order of magnitude in excess of the theoretical prediction based on the CSM. 

In order to convincingly establish the phenomenological signihcance of the CO pro¬ 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and vN deep-inelastic 
scattering, annihilation, 77 collisions, and 6 -hadron decays may be found in the 

literature; see Ref. and references cited therein. Furthermore, the polarization of char¬ 
monium, which also provides a sensitive probe of CO processes, was carefully investigated 
P JT0|JTl| . None of these studies was able to prove or disprove the NRQCD factorization 
hypothesis. 

Very recently, the DELPHI Collaboration has presented preliminary data on the in¬ 
clusive cross section of J/ijj photoproduction in 77 collisions —*• e^e~ J/p + X) at 

CERN LEP2, taken as a function of the J/p: transverse momentum [px) ig. The Jli< 
mesons were identihed through their decays to pairs, and events where the system X 
contains a prompt photon were suppressed by requiring that at least four charged tracks 
were reconstructed. The luminosity-weighted average e’''e“ center-of-mass (c.m.) energy 
was ^/s = 197 GeV, the scattered positrons and electrons were antitagged, with maximum 
angle ^max = 32 mrad, and the 77 c.m. energy was constrained to be IT < 35 GeV in 
order to reject the major part of the non-two-photon events. The total cross section was 
found to be cr(e+e“ — »• e^e~ J/ip + X) = (45.3 ± 18.8) pb. 

In this Letter, we seize the opportunity to confront this data with up-to-date theo¬ 
retical predictions based on NRQGD and the GSM, in order to hnd out if it is able to 
discriminate between the two. In want of the full next-to-leading-order (NLO) corrections, 
it is indispensable to perform a comprehensive and conservative analysis of the theoret- 
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ical uncertainties. Doing this, we shall hnd that this data clearly favors the NRQCD 
prediction, while the CSM one signihcantly falls short of it. 

The photons can interact either directly with the quarks participating in the hard- 
scattering process (direct photoproduction) or via their quark and gluon content (resolved 
photoproduction). Thus, the process e’''e“ —*• e^e~J/'ilj + X receives contributions from 
the direct, single-resolved, and double-resolved channels. All three contributions are for¬ 
mally of the same order in the perturbative expansion and must be included. This may 
be understood by observing that the parton density fnnctions (PDFs) of the photon have 
a leading behavior proportional to oc a/as, where a is Sommerfeld’s 

hne-strnctnre constant, M is the factorization scale, and Aqcd is the asymptotic scale 
parameter of QCD. 

In 77 collisions, J/'ip mesons can be produced directly; or via radiative or hadronic 
decays of heavier charmonia, snch as Xcj and X mesons; or via weak decays of b hadrons. 
The respective decay branching fractions are B{xco = (0.66±0.18)%, B{xci 

J/^P + X) = (27.3 ±1.6)%, B{xc2 ^ W + 7) = (13.5 ±1.1)%, S(V’' ^ J/V’±X) = (55 ± 
5)%, and B{B —>■ J/'ip + X) = (1.16 ±0.10)% [0. The OPAL Collaboration has recently 
measured the total cross section of open-bottom production in 77 collisions at LEP2, 
nnder similar kinematic conditions as DELPHI [|T^] (a/s = 194 GeV, = 33 mrad, 
and 10 < hF < 60 GeV), to be a{e^e~ —>• e~^e~bb + X) = (14.2 ± 5.9) pb |]T4[. The cross 
section for J/'ip mesons from 6-hadron decays may thus be estimated to be (0.33±0.14) pb, 
which is less than 1% of the total J/ip cross section measured by DELPHI and can be 
safely neglected. The cross sections of the fonr residnal indirect prodnction channels may 
be approximated by mnltiplying the direct-prodnction cross sections of the respective 
intermediate charmonia with their decay branching fractions io J/ip mesons. 

Invoking the Weizsacker-Williams approximation 0 and the factorization theorems 
of the QGD parton model and NRQGD [Q], the differential cross section of e+e“ —> 
e^e~H ± A, where H denotes a generic charmonium state, can be written as 


da{e^e —^ e^e H + X) = J dx+fj/e{x+) j dx- 

X f^/e{X-) / dXafa/-yiXa, M) / dXbfb/j{Xb, M) 

a,b,d'' '' 

[n])da{ab —>• cc[n] ± d), 


( 1 ) 


where f^/eix±) is the eqnivalent number of transverse photons radiated by the initial- 
state positrons and electrons 0 . fa/'y{xa, M) are the PDFs of the photon, {0^[n\) are 
the MEs of the H meson, da{ab —>■ cc[n] ± d) are the differential partonic cross sections, 
the integrals are over the longitudinal-momentnm fractions of the emitted particles w.r.t. 
the emitting ones, and it is snmmed over a, 6 = 7 , g, q, q and d = g,q, q, with q = u,d, s. 
To leading order in v, we need to include the cc Fock states n = 
if 77 = J/ip,ip' and n = if 77 = XcJi where J = 0,1,2. With the dehnition 

f^/^{xy, M) = 5(1 — Xy), Eq. dH) accommodates the direct, single-resolved, and donble- 
resolved channels. The presence of parton d is to ensure that pt can take hnite valnes; 
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if it were absent, then px would essentially be zero, so that only the lowest bin of the 
DELPHI data could be described. 

We analytically calculated the cross sections of all contributing partonic subprocesses 
and compared our results with the literature.[] Specihcally, these include 77 —> cc[^S[^'^]g 




cc[^Pj^^^ 


cc[^S[% [|, cc[8]g 
cc[8]g §, 0 ; gq 


7 g —>• cc[ 8 ]q' 
cc[8]q 


21 


99 


cc[^S[^^]g 
>(i)i 


’]9 1, cc[ 8 ]^ §,0; gq cc[^Py]q [|, cc[ 8 ]g [|,|T0l; and qq [^Pr]g 
cc[8]g [|6|,|l0i, where n = 8 collectively denotes n = ^S^q\^s[^\^Pj‘\ In the limit px —> 0, 
some of these cross sections are plagued by collinear and infrared singularities. This 
happens whenever the respective 2 —>• 1 partonic subprocess, ab —>• cc[n], exists, namely, 
for 77 —>• ccpPo 2 ^]) 79 cc[ 8 '] gg —>• cc[8'] and qq —>• ceps'®] where 


n = 8' stands for n = ^S*!] ,^PoJ- In a full NLO analysis, these singularities would be 
factorized at scale M and absorbed into the bare PDFs and MEs so as to renormalize the 
latter. In our LO analysis, we refrain from presenting predictions for the px distribution in 
the lowest px bin. Instead, we consider the cross section arising from the 2 —> 1 partonic 
subprocesses. 

In our numerical analysis, we use me = (1.5 ± 0.1) GeV, a = 1/137.036, and the 
lowest-order (LO) formula for as^^\p,) with Uf = 3 active quark flavors. As for the 
photon PDFs, we use the LO set from Gliick, Reya, and Schienbein (GRS) ^ 5 , which 
is the only available one that is implemented in the hxed-flavor-number scheme, with 
Uf = 3. We choose the renormalization and factorization scales to be /r = ^mx and 

M = respectively, where mx = \J{2my + p\ is the transverse mass of the H 

meson, and independently vary the scale parameters y and between 1/2 and 2 about 
the default value 1. As for the J/ip, Xcj, and y ME’s, we adopt the set determined 
in Ref. by htting the Tevatron data using the LO proton PDFs from Martin, 
Roberts, Stirling, and Thorne (MRST98LO) as our default and the one referring to 
the LO proton PDFs from the GTEQ Gollaboration (GTEQ5L) for comparison (see 
Table I in Ref. [0). In the hrst (second) case, we employ Aq^d = 204 MeV (224 MeV), 
which corresponds to Aq^^ = 174 MeV (192 MeV [^), so as to conform with the 
ht |]TI|. Incidentally, the GRS photon PDFs are also implemented with Aq^^ = 204 MeV 
In the cases -i/ = , the ht results for ((P’^[^S'®]) and {0^[^Pq^'^]) are strongly 


correlated, and one is only sensitive to the linear combination 


Mp = (O^pA®]) + (2) 

with an appropriate value of r [|6|, pR| , pd| . Since Eq. (|l|) is sensitive to a different linear 
combination of (C>^[^S'®]) and than appears in Eq. (||), we write ((P’^[^S'®]) = 

nMf and ((P’^pPo^^j) = (1 — k) (ml/r) Mf and vary k between 0 and 1 about the default 
value 1/2. The J-dependent MEs ((P’^pP®]), ((P^'=-^pPyp), and ((P^'^-^pS'®]) satisfy the 
multiplicity relations collected in Eq. (4) of Ref. |^, which follow to leading order in v 

^Leaving aside obvious typographical errors, we disagree with Eq. (7) and the subsequent equation in 
Ref. @ and with Eq. (4) in Ref. ||^. 
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from heavy-quark spin symmetry. In order to estimate the theoretical uncertainties in 
our predictions, we vary the unphysical parameters and k as indicated above, 

take into account the experimental errors on rric, the decay branching fractions, and the 
default MEs, and switch from our default ME set to the CTEQ5L one, properly adjusting 
Aqqd- We then combine the individual shifts in quadrature, allowing for the upper and 
lower half-errors to be different. 

In Fig. m, we confront the distribution of —*• e^e~ J/'ip + X measured by 

DELPHI with our NRQCD and CSM predictions. The solid lines and shaded bands 
represent the central results, evaluated with our default settings, and their uncertainties, 
respectively. We observe that the DELPHI data clearly favors the NRQCD prediction, 
while it signihcantly overshoots the CSM one. This is even more apparent from the 
data-over-theory representations shown in Fig. This qualitative observation can be 
substantiated by considering the values for the N = 9 data points with > 0.25 GeV^. 
In fact, the NRQCD central prediction yields x^/iV = 0.49, which is to be compared 
with 1.79 for the CSM one. The situation is very similar for the MEs pertinent to the 
CTEQ5L PDFs (dashed lines), the corresponding results being 0.62 and 1.76, respectively. 
As for the integral over the range 1 < < 10 GeV^, the DELPHI, NRQGD, and 

GSM results read (6.4 ± 2.0) pb, pb, and 0.39(“jQQg pb, respectively, he., the 

DELPHI measurement and the NRQGD prediction mutually agree within errors, while 
the GSM prediction signihcantly falls short of the DELPHI result, by a factor of 16 as 
far as the central values are concerned. In Fig. ^ we study the normalized cross section 
{I/a)da/dp"^, which is particularly sensitive to the shape of the distribution and offers 
the advantage that the theoretical uncertainty is greatly reduced. We observe that both 
NRQGD and the GSM describe the shape of the measured p^ distribution well within its 
errors, the respective x^/-^ values being 0.35 and 0.68, respectively. 

Taking as the reference quantity the p\ distribution integrated from 0.25 to 10 GeV^, 
the direct, single-resolved, and double-resolved channels account for 1%, 98%, and 1% of 
the NRQGD prediction, respectively. The situation is very similar for the GSM, except 
that direct photoproduction is forbidden at LO. In NRQGD, 91% of the J/'ip mesons are 
directly produced, while 9% stem from the decays of Xcj and -0^ mesons. In the GSM, 
direct production happens less frequently, in 77% of the cases. As explained above, the 
contribution from 6 -hadron decays is negligible. Gonsequently, the ratio of the direct and 
indirect yields lends itself as a useful discriminator between NRQGD and the GSM, and 
it would be desirable to measure it in 77 collisions. In NRQGD, the most important error 
sources include the variations of rUc, and k, which, in average, make up 55%, 24%, 
and 11% of the total error square + ( 6 |Ct)^, respectively. In the GSM, the largest 

errors are related to (59%), tjIc (30%), and P,m (7%). It is generally believed that the 
magnitude of unknown higher-order corrections may be estimated by scale variations of 
the known results. For the single-resolved channel, which greatly dominates our results, 
the NLO corrections in the GSM are known to be moderate at small values of pt |2^. 
A similar observation was made in Ref. I12D 


for the direct channel, by studying the real 
radiative corrections in full NRQGD. This indicates that the and variations of the 
LO result may indeed provide a realistic estimate of the size of unknown higher-order 
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corrections. 

Recent analyses indicate that the overall description of the Tevatron |0 and HERA 
data on inclnsive charmoninm prodnction can be improved by adopting the kx factoriza¬ 
tion approach. It wonld be interesting to hnd ont if this is also trne for the DELPHI data 

0 - 
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Figure 1: The cross section da/dp^ of e'''e“ —>■ e^e~ J/ip + X measured by DELPHI [H 
as a function of is compared with the theoretical predictions of NRQCD and the CSM. 
The solid and dashed lines represent the central predictions obtained with the ME sets 
referring to the MRST98LO (default) and CTEQ5L PDFs, respectively, while 
the shaded bands indicate the theoretical uncertainties on the default predictions. The 
arrows indicate the NRQCD prediction for pr = 0 and its ^S[^\ and 

components. 
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Figure 2: Data-over-theory representation of Fig. |l|. 
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Figure 3: The normalized cross section {l/a)da/dp^ of e'''e“ —»• e~^e~J/ip + X measured 
by DELPHI as a function of is compared with the default theoretical predictions 
of NRQCD and the CSM. The shaded bands indicate the theoretical uncertainties. 
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